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1. Introduction 

Let us consider the Green function for the steady-state Schrodinger equation in one 
dimension, 

--^^i^) + VsixMx) = k^i^ix). (1.1) 

The Green function plays important roles in various physical problems, and there are 
many approaches to the study of the Green function. In this paper we discuss a new 
description of the Green function in terms of reflection coefficients. 

From the physical point of view, it is natural to interpret the propagation of waves in 
terms of the processes of multiple reflections and transmissions. In quantum mechanics, 
this interpretation has been used, for the most part, in the context of semiclassical 
approximations [1, 2]. The Bremmer series, which is a perturbative improvement 
of the WKB approximation, is based on this picture [3]. A similar idea is used in 
the invariant imbedding method [4], which has applications in many areas including 
transport problems in astrophysics, conductors, and random media [5-8]. The essence 
of the invariant imbedding method is to express everything in terms of "emergent" or 
"observable" quantities such as transmission and reflection coefficients, without need 
of considering what is happening within the system. In this method one deals with 
reflection coefficients for finite intervals, and, by varying the endpoint of the interval, 
derives a differential equation of Riccati type satisfied by the rcfiection coefficients. The 
derivation of this Riccati equation is essentially equivalent to taking account of the 
transmission and reflection processes at the endpoint. 

It is possible to use the same idea to construct the Green function. By taking 
the sum over all the multiple reflections and transmissions, we can derive exact 
expressions for the Green function [9, 10] J. These expressions are written in terms of 
the transmission coefficient for a finite interval, and the reflection coefficients for finite 
and semi-infinite intervals. With these expressions, the analysis of the Green function 
can be reduced to that of the transmission and reflection coefficients. 

The structure of reflection coefficients for semi-infinite intervals have been throughly 
studied, and various formulas have been obtained for their high- and low-energy 
behaviors [11]. However, the mathematical structure of transmission coefficients is 
not as simple. This is because transmission coefficients are "non-local" quantities 
in the sense that they are functions of two endpoints of the finite interval. We 
may say that, in a sense, reflection coefficients are more fundamental quantities than 
transmission coefficients. The analysis of the Green function becomes much easier if it 
is expressed solely in terms of reflection coefficients for semi-infintc intervals, without 
using transmission coefficients. It is the objective of the present paper to derive such 
expressions. 

I The expression for the Green function derived in [10] for segmented potentials is identical to the one 
obtained in [9] for the Fokker-Planck equation. 
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The expressions in terms of reflection coefficients are particularly useful for the 
analysis in the high- and low-energy regions. By using the formulas already known for 
the reflection coefficients, we can derive new formulas for asymptotic expansions of the 
Green function. The advantage of this approach over conventional methods is that it 
can be applied to a larger class of potentials. The reflection coefficients can be defined 
irrespective of whether the potential Vs{x) is finite or infinite as x — > ±00; we do not 
need to assume that Vs (x) vanishes sufficiently rapidly at infinity, as is necessary for the 
description using Jost solutions. We do not need to care about the existence of bound 
states, nor do we have to know the eigenvalues. Conventional methods of analysis 
are sensitive to the behavior of the potential at infinity, and it is often necessary to 
use different methods for different kinds of potentials. In the formulation in terms of 
reflection coefficients, the analysis of the Green function can be carried out for various 
types of potentials in a unified way. In addition, the formulas for asymptotic expansions 
obtained in this method are more explicit than the ones obtained by conventional 
methods. (This will be discussed in a separate paper.) 

The expressions in terms of reflection coefficients are also convenient for calculating 
the Green function in practical situations, either approximately or numerically. It 
turns out that the expressions derived in this paper have a close relation with the 
WKB method. In the light of the formalism developed here, we can understand the 
WKB method from a new viewpoint, which may possibly lead to new improvements 
of the WKB approximation. Our expressions can also be used as a basis for other 
new approximation methods. Since the reflection coefficients are quantities that have 
a clear physical meaning, expressing the Green function in terms of them is useful 
for the purpose of making approximations. The reflection coefficients are also suited 
for numerical treatments, and so these expressions will be useful for the numerical 
calculation of the Green function, too. 

In our method of derivation, we make use of the Fokker-Planck equation. It is 
well known that the Schrodinger equation (1.1), with an appropriate shift of the energy 
level, can be transformed into a Fokker-Planck equation [12]. The (time-independent) 
Fokker-Planck equation describing the Brownian motion in a potential V{x) has the 
form 

--^m + 2^[f{xmx)] = k'cpix), (1.2) 

where 

Equation (1.2) is equivalent to (1.1), where 

^^;{x) = e^(-)/20(x), Vs{x) = f{x) + f{x). (1.4) 

With the use of the Fokker-Planck equation, it becomes easier to study the structure 
of the transmission and reflection coefficients, and various formulas take simpler forms. 
In particular, a symmetry transformation of the Fokker-Planck equation plays a crucial 
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role in our method. As a result we obtain a onc-paramcter family of expressions, which 
reflects the symmetry structure of the Fokker-Planck equation. 

We assume that Vs(x) either converges to a finite value or diverges to +00 as 
X — )■ +00, and that Vs(a;) is also either finite or +00 as a; — )■ —00. (We do not consider 
the cases where Vs{x) tends to —00 as x ^ ±00, or the cases where Vs{x) oscillates at 
infinity.) We also assume that k is, in general, a complex number with ImA; > 0. Let 
Gs{x,x';k) denote the Green function for equation (1.1), satisfying 

Gs{x,x';k)^5{x-x') (1.5) 

with the boundary condition Gs{x, x'\ /c) ^ as |a; — x'\ — >■ 00 for Im A; > 0. We define 

G{x, x'; k) = 2ikGs{x, x'; k). (1.6) 

In this paper we shall deal with the quantity G defined by (1.6), rather than Gs itself. 
(For convenience, we shall also call G the Green function.) Without loss of generality 
we may assume that x > x'. The expressions for x < x' are obtained by interchanging 
X and x'. 

Let us now define the refiection coefficients for semi-infinite intervals. For general 
Vs{x) (rather than special forms such as piecewise constant or segmented potentials), 
there is no unique natural way of defining the refiection coefficients for finite or semi- 
infinite intervals. As mentioned above, we shall define them in terms of the Fokker- 
Planck equation, and this turns out to give the simplest description. Our definition of 
the reflection coefficients for semi-inflntc intervals is illustrated in flgurc 1. Let Xq be an 
arbitrarily chosen point. We let the Fokker-Planck potential V{x) in the region x > Xq 
be replaced by the constant value V{xq), and dcflnc 

V{x) = V{x)e{xo -x) + V{xo)e{x - xo), (1.7) 

where 9 is the Heaviside step function. (Recall that the Schrodinger potential Vs{x) is 
related to V{x) by equations (1.3) and (1.4).) We consider equation (1.2) with f{x) 
replaced by f{x) = —{l/2){d/dx)V{x). In the region x > Xq, where f{x) = 0, this 
equation has independent solutions of the form e^*'^^ and e~**^^. We deflne the reflection 
coefficient Rr{xo, —oo;k) as the coefficient multiplying the reflected wave e**^(^~^o) in the 
region x > xq when there is an incident wave e"*'^^^"'^''). In other words, Rr{xo, —00; k) 
is defined by a solution of the form 

(f){x) = e-'^^^^-^o) + R,{xo, -00; A;)e'*=(^-^°) for x > xq, (1.8a) 
0(x) ^0 as X ^ -00. (1.86) 

(When k is real, it is necessary to assume in (1.86) that k has an infinitesimal imaginary 
part ie with e > 0.) If Vs{x) = for < a; with some Xi, and if Xi < Xq. then the 
above definition of Rr{xo, —00; k) coincides with the usual definition of the reflection 
coefficient. In general cases, the Schrodinger potential corresponding to the Fokker- 
Planck potential (1.7) includes a delta function ect x — xq- In the same way, the left 



dx' 



-Vs(x) + k' 



Expressions of the Green function 



5 



V{x) 




Figure 1. Definition of Rr{xo, — oo; k) and Ri{oo, xq; k). 

reflection coefiicient for the interval {xq, +oo) is deflned by considering, instead of (1.7) 
and (1.8), 

V{x) = V{x)0{x - xo) + V{xo)e{xo - x), (1.9) 

and 

(f){x) = e^*=(^-^o) + Ri{oo, Xo; A;)e-^'=(™) for x < xo, (1.10a) 
(t){x) as X +00. (1.106) 

Our objective is to express the Green function in terms of these two quantities, 
Rr{xo, — oo; A:) and Ri{oo, xq; k). The results arc applicable to any Vs(a;) (which is either 
finite or +oo at a; = ±oo) as long as the refiection coefficients can be defined for it. 

2. Boson representation 

It was shown in [13] that the Green function can be expressed in a general form in 
terms of the Lie superalgebra osp{l/2). We can obtain various expressions of the Green 
function by writing this general expression in specific representations. Here we use a 
representation in terms of boson operators, which is convenient for the methods we shall 
use in this paper. 
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Let a and be the boson annihilation-creation operators, satisfying the 
commutation relation 

[a,at] = l, (2.1) 

and let |0) be the boson vacuum state, satisfying 

a|0)=0, (0|at = 0, (0|0) = 1. (2.2) 

We regard the space coordinate x as playing the role of the time, and consider the 
"Hamiltonian" 

H(x) = -k (a^a + |) - ^if(x)(aa - a^a^), (2.3) 

where f{x) is the function defined by (1.3). The free part of this hamiltonian, —ka^a 
describes free propagation of the boson. The interaction part consists of pair creation 
and pair annihilation of bosons. (The constant term —^k is added for later convenience.) 
We define the evolution operator U (x, xq) as the solution of the differential equation 

d 

i—U{x, Xq) = H{x)U{x, xo) (2.4) 
ox 

with the initial condition U{xo,Xo) = 1. Using this evolution operator, G{x,x') defined 
by (1.6) can be written as [14] 

, {0\U{^,x){a + a^)U{x,x'){a + a^)U{x',-^)\0) 

= (0|^(oo,-oo)|0) ■ ^'-'^ 

This is a specific form of the general algebraic expression mentioned above. To 
understand the meaning of this expression, it is helpful to think about the expansion 
of the right-hand side in powers of /. This expansion can be visualized by using 
Feynman diagrams. Graphically, the right-hand side of (2.5) is obtained as the sum 
of all connected diagrams like the one shown in figure 2(a). (Disconnected diagrams are 
cancelled by the vacuum amplitude in the denominator.) Each diagram represents a 
path connecting the points x' and x. The rules for interpreting the diagrams are shown 
in figure 2(b). It should be noted that the expression (2.5) itself is valid even when the 
expansion in terms of / is not well defined, e.g., when f{x) is infinite ai x — ±oo. 

3. Scattering coefficients and the Green function 

The scattering coefficients for a finite interval {xi,X2) arc defined in the same way as the 
reflection coefficients for semi-imfinte intervals we have already introduced. We consider 
the Fokker-Planck potential 

{V{xi) X < Xi, 
V{X) Xi<X< X2, (3.1) 

V{X2) X2 < X. 

Equation (1.2) with f{x) replaced by f{x) = —{l/2){d/dx)V{x) has two independent 
solutions of the form 

. e[^(^^)-^(^i)l/V(x2, XV, A;)e-^'=(^-^i) x < x,, 

M^) ^ e-ik{x-X2) ^j^^^^^^^^.}^yk{x-x,) X2<X, 




g *(Z2-Zl ) 

(b) ; , 

Zi Zi 

Figure 2. (a) A typical diagram connecting the points x' and x. (The vertical 
direction of this figure does not have any particular meaning.) (b) The diagrammatic 
rules. Each line segment connecting xi and X2 corresponds to the free propagator 
^ik(x2-xi) ^ To each turning point of the path is assigned a factor ±/(z), where the 
sign is plus if the path comes to that point from the right, and minus if it comes from 
the left. 




e'^(^-^^) + Ri{x2. X,- fc)e-^'=(^-^^) x < x,, 

M^) - S -[y(..)-y(..)]/2^(^^^ fc)e^^(— < x. ^ ^ ^ 

This defines the transmission coefficient r, the right reflection coefficient Rr, and the 
left reflection coefficient Ri for the interval (xi,X2). In the boson representation, they 
can be written as [14] 

{0\aU{x2,Xi)a^0) 

"'^'■"''^ (o|c/(.....)|o> ■ 

(0|f/(X2,Xi)|0) 

(0|f/(x2,xi)atat|0) 
(0|^(x2,x,)|0) • ^'-''^ 

The expressions (3.36) and (3.3c) also hold for semi-inflnte intervals. Namely, 
Rr{xo, — oo) and -R;(+oo,xo) defined by (1.8a) and (1.10a) can be expressed as 

(O|aa[/(xo,-oo)|0) (0|^(oo, Xo)atat|o) 

Rr{xo, -oo) = --— , Ri{oo, xo) = (3-^) 

(0|f/(xo, -oo)|0) (0|[/(cx),Xo)|0) 

Similarly to the graphical interpretation of G{x, x') shown in figure 2, we can interpret 
(3.3) in terms of diagrams. The transmission coefficient t{x2,Xi) is the sum of all 
the paths that penetrate the interval (xi,X2), as in figure 3(a). The diagrams for the 





Xl X2 



Figure 3. A typical diagram or (a) r, (b) Rr, and (c) Ri. Sucli diagrams are to be 
evaluated with the rules given in figure 2(b). 



reflection coefficients consist of tlie patlis tliat start from one of the endpoints of the 
interval and return to that same point, as shown in figures 3(b) and 3(c). 

The scattering coefficients for finite intervals are the quantities that play major 
roles in the invariant imbedding method. We shall use them as building blocks for 
constructing the full propagator (2.5). However, these quantities shall appear only 
in intermediate steps and not remain in our final results. Our objective is to express 
everything in terms of the reflection coefficients for semi-infinite intervals, without using 
the quantities (3.3) for finite intervals. 

As explained in section 2, the propagator G{x, x') is the sum of the paths connecting 
the points x' and x. Such paths can be constructed from the transmission and reflection 
coefficients. The idea used here is essentially the same as the old one which dates back 
to the work by Stokes [15]. As illustrated in figure 4, we have [9] 

G{x,x') = [1 + Ri{oo,x)](f2[Ri{oo,x)Rrix,-oo)]A 

X r(x, x')(Y.[Ri{x, x')Rrix', -oo)]'"] [1 + i?,(x', -oo)]. (3.5a) 

Note that this expression is not symmetric with respect to x and x'. This asymmetrical 
treatment is necessary in order to avoid double counting. It is also possible to exchange 
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Figure 4. The passage from x' to x can be decomposed into the processes shown 
here. They correspond to: (a.) 1 + Rr{x' , —oo), (b) ^^Q[i?;(a;, a;')i?r(a;', — oo)]", (c) 
t{x,x'), {d) Y,m=oi^i-i°°^^)^r{x,-oo)]"', {e) 1 + Ri{oo , x) . 



the roles of x and x' in (3.5a) and write 

^ oo 

G{x, x') = [1 + Ri{oo, x)] ( ^ [Riioo, x)Rr{x, x')y 



X t{x, x') ( £ [Riioo, x')Rr{x', -oo)]"''] [1 + Rr{x' , -oo)]. (3.56) 

\n=0 ' 

The geometric series in equations (3.5) can be summed to yield 

[1 + i?i(oo,x)][l + i?^(x', -oo)]r(x,x') 

G(X,X) = -f ; ; — ; ^- ; — , (3.6a) 

[1 — Ri{oo, x)Rr{x, — oo)J[l — Ri{x, x')Rr{x', —oo)\ 

, ^ [1 + Rijoo, x)] [1 + Rrjx', -oo)]r(x, x') 
^^'^ ^ ~ [1- Ri{oo,x)Rr{x,x')][l- Ri{oo,x')Rr{x',-oo)]' ^' ' 

We wish to eliminate the r(x, x'), Rr{x,x'), and Ri{x,x') from these expressions. We 
shall do this in the next section. 



4. Expressions in terms of reflection coefficients 

The vacuum amplitude {0\U{x2, xi)|0), which appears in the denominators on the right- 
hand sides of (3.3), is related to the transmission coefficient by the identity [13] 

{0\U{x2,x,)\0) = [r{x2,x^)f\ (4.1) 

This identity can be checked diagrammatically for each order in / (figure 5). Since 
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Figure 5. Diagrammatic interpretation of the identity (4.1) to (a) order and (b) 
order On the left-hand sides are the the diagrams of r. On the right-hand sides, 
the straight hnes correspond to e'^'^'^i^^-^i) ^ which comes from the constant term —-^k 
in the Hamihonian. The bubbles are the diagrams of Z^, where Z is defined by (4.2). 

(Using the quantity F defined by (4.3), we can write Z^ = 1 + 2F + {2Fy/2 H . As 

shown in figure 6, the diagrams for F consist of connected loop diagrams.) 



a constant term —^k is included in the Hamiltonian (2.3), we have (0|?7(x2, Xi)|0) = 
gifc(a;2-xi)/2 -^j-^gj^ J jg identically zero. We define 

Z{x2,x,) = e-^^(^2-^'^)/2^0|t/(x2,a;i)|0), (4.2) 

so that Z = 1 when / = 0. This Z is the vacuum amplitude in the usual sense; the 
Feynman diagrams for Z are bubble diagrams without external legs. These bubble 
diagrams are, in general, disconnected. To deal with connected Feynman diagrams, we 
define 

F(x2, xi) = log Z{x2, xi). (4.3) 

As is known in usual diagrammatic discussions in field theory [16], this F is obtained as 
the sum of all connected loop diagrams (figure 6(a)). (In statistical mechanics, Z and 
F correspond to the partition function and the free energy, respectively.) From (4.1), 
(4.2), and (4.3), we have 

r(x2, xi) = exp [ik{x2 - Xi) + 2F{x2, Xi)] . (4.4) 

Connected loop diagrams are obtained by connecting the two legs of Rr with a 
factor — / (see figure 6(b)). This fact can be expressed as 

1 /■^2 

F{x2,xi) = -- f{z)Rr{z,xi) dz. (4.5a) 

Z Jxi 

(There is a factor 1/2 because the same diagram is obtained by exchanging the two legs 
oi Rr.) In the same way, F can also be obtained from Ri as shown in figure 6(c). So we 
have 

1 /"^a 

Fix2, xi) = - f{z)Ri{x2, z) dz. (4.56) 




Figure 6. (a) A typical diagram of F{x2,xi). As shown in (b) and (c), such a diagram 
can be obtained from a diagram for or Ri (figures 3(b) and (c)). 



In the invariant imbedding method, differential equations satisfied by the scattering 
coefficients are derived by varying an endpoint of the interval. Analogous differential 
equations for the quantity F are obtained from equations (4.5) as 

d 1 d 1 

F{x2,Xi) = --f{x2)Rr{x2,Xi), -— F {X2, Xi) = - - f (xi) Ri{x2, Xi) . (4.6) 



dx2 2''-' dxi 2 

There is another useful relation that connects F to the reflection coefficients: 

F{Xc, Xa) - F{Xc, Xb) - F{Xb, Xa) = log [1 - Rl{Xc, Xb)Rr{Xb, Xa)] , (4.7) 

where < x;, < Xc- We can understand this relation diagrammatically. The left-hand 
side of (4.7) is the sum of all connected loop diagrams which are restricted within the 
interval {xa,Xc), and which cross the point Xb (figure 7(a)). Such diagrams are obtained 
from the reflection coefficients as shown in figure 7(b). The series in figure 7(b) can be 
summed as 

RlRr + \{RlRrf + \{RlRrf + ■■■ = - log(l - RiRr). (4.8) 

Hence we have (4.7). (There is an overall factor ^ on the right-hand side of (4.7) for 
the same reason as in equations (4.5).) All the relations such as (4.5) or (4.7) are valid 
even when the expansion in terms of / is not well defined. (It is not difficult to prove 
these relations without using the diagrams.) 

Differentiating both sides of (4.7) with respect to Xa, Xf,, or Xc, and using (4.6), we 
obtain 
d 

— log [1 - i?z(Xc, Xb)Rr{Xb, Xa)] = f{Xc) [Rr{Xc, Xa) - Rr{Xc, Xb)] , (4.9a) 

d 



dx} 



log [1 - Rl{Xc,Xb)Rr{Xb,Xa)] = - f {Xb) [Rl{Xc,Xb) + Rr{Xb,Xa)] , (4.96) 
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(b) 



R, 



+ 



R, 



R, 









/ 











Ri 



Ri 



+ I 



Rr 




R, 




Rr 


^' 


Ri 


+ 


Rr 




Ri 





Figure 7. (a) A diagram contributing to F{xc, Xa) — F{xc, Xb) — F{xi,, Xa). (b) Such 
diagrams can be constructed in this way, using the reflection coefficients. (The factors 



2' 3 



etc are necessary in order to avoid double counting.) 



d 

dXa 



log [1 - Rl{Xc, Xb)RriXb, Xa)] = f{Xa) [Rl{Xc, Xa) - Rl{Xb, Xa)] ■ 



(4.9c) 



Setting Xc = z, Xb = xi, Xa = — oo in (4.9a), and integrating both sides with respect to 
z from Xi to X2, we have 

\og[l - Ri{x2,xi)Rr{xi,-oo)] = f{z)[Rr{z,-oc)-Rr{z,xi)]dz, (4-10) 

J Xl 

where we have used Rr{xi,Xi) = 0. Using (4.10), we can rewrite (4.5a) as 

F{x2,xi) = -- f{z)Rr{z, -00) dz + -\og[l- Ri{x2,Xi)Rr{xi, -00)]. (4.11) 
Z Jxi z 

Substituting this into (4.4) yields 

t{x^x') = [1 — Ri{x,x')Rr{x',— 00)] exp ik{x — x') — / f{z)Rr{z, —00) dz . (4.12a) 

L Jx' J 

In the same way, using (4.56) and (4.9 c) we can derive 



T[X, X 



[1 — Ri{oo, x)Rr{x, x')] exp 



ik{x — x') + / f{z)Ri{oo,z)dz 



(4.126) 



On substituting (4.12a) into (3.6a), the factor including Ri{x,x') cancels out, and 
G{x, x') is expressed solely in terms of reflection coefficients for semi-infinite intervals: 



[1 + Ri{oo,x)][l + Rr{x',-oo)] 

G{x, X = — — ^ exp 

1 — Ri[oo, x)Rr{x, —00) 



ik{x — x') — I f{z)Rr{z,—(yo)dz . 

(4.13a) 
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Similarly, from (4.12&) and (3.66) we obtain 



G(x, x') = 



[1 + Ri{oo, x)][l + Rrjx', -oo)] 
1 — Ri{oo, x')Rr{x', — oo) 



exp 



ik{x — x') + I f{z)Ri{oo,z)dz 

Jx' 



On the other hand, integrating (4.96) and setting Xc — oo, Xa — — oo gives 

""fl WUt A, D/ 1 I- Rl{^,x')Rr{x',-(X>) 

f[z)[Ri[oo, z) + Rr[z,—oo)\ dz = \og- — 

x' 



(4.136) 



(4.14) 



' ^ ' 1 — Ri{oo,x)Rr{x,—oo) 

It is obvious that (4.13a), (4.136), and (4.14) are consistent. A symmetric expression of 
G{x,x') is obtained by multiplying (4.13a) and (4.136), and taking the square root: 

[1 + i?i(oo, a;)][l + i?^(a;', -oo)]e^*^(^-^') 



G{x,x') 



[1 - Ri{00, x)Rr{x, -oo)] V2[l _ i?;(oo, x')Rr{x', -Oo)]V2 



(4.15) 



5. Generalization 



In this section we shall derive a more general expression of the Green function, which 
includes (4.15) as a special case. The derivation is based on a symmetry transformation 
which can be understood as a rotation of the coordinate axes [17]. 
We define 

X{x) = ikx, Y{x) = V{x)/2, (5.1) 
where V{x) is the Fokker-Planck potential. Then (2.4) with (2.3) can be written as 

g^U{x,Xo) = - 

We consider the rotation of the X-Y axes, defining 

f Xe\ f cose sine \ ( X \ 
[ Yn I ~\ -sine cose I [ Y j' 



dX i J. dY t t\ 

— — (aa' + a'a) H — ;— (aa — a'a'j 
dx dx 



U{x,xo). 



(5.2) 



(5.3) 



We also define the boson operators in the rotated frame as 



J 



cos I — sm 2 
sin I cos I 



a 
at 



(5.4) 



They are indeed boson operators, satisfying the commutation relation 

[ae,aj] = l. (5.5) 

It is easy to see that equation (5.2) is covariant under this rotation; it sill holds when 
X, Y, a, are replaced by the ones with subscript e: 



^Uix,x,)^l 



Let |0; e) denote the vacuum state in the rotated frame, satisfying 
a(^|0; e) = 0, (0; e\al = 0, (0; ^|0; e) = 1. 



U{x,Xo). (5.6) 



(5.7) 
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Z2 



Figure 8. The diagrammatic rules in the rotated frame with angle 9. 



It can be shown that this state is related to the original vacuum as [14] 

\0-e) = {l + 7]^y/^exp{7]a^ay2)\0), (5.8a) 

(0; ^1 = (1 + r/2)^/^(0| exp{-ri aa/2), (5.86) 

|0) = (1+7^2)1/4 exp(-r^44/2)|0;^), (5.9a) 

(0| = (1 + 7/2)1/4(0; ^|exp(r7aeae/2), (5.96) 

where 

r/ = tanf. (5.10) 

Using (5.9), we can rewrite (2.5) as 

_ (0;e|e''"«'^''/2f/(oo,a;)(a + at)?7(x,x')(a + at)f/(a;',-cx))e~''"9"«/2|o.^) 



(0; ^|e''»«"''/2?7(oo, -oo)e-i''l''e/^\0; 9) 



(5.11) 



where 

a + = {cos I — sin |^ ae + (cos | + sin |^ a\. (5-12) 

Equation (5.11) holds for any 9, and so it is a generalized form of (2.5). 

Just like (3.3), we define the scattering coefficients in the rotated frame: 

{Q]9\aeU{x2,Xi)al\Q-9) /c- io n 

reX2,Xi = — — — — , 5.13a 

(0; 9\U[x2,xi)\Q] 9) 

{Q]9\aeagU{x2,Xi)p,9) /crion 

, . _ (0;g|[/(3:2,3;i)a^at|0;g) 
(0;6'|t/(x2,xi)|0;6') 

Since the evolution equation (5.6) has the same form as (5.2), these scattering coefficients 
can be interpreted diagrammatically in the same way as before (i.e., as in figure 3). The 
rules in figure 2(b) are now generalized to the ones shown in figure 8; as can be seen 
from the right-hand side of (5.6), the free propagator connecting xi and X2 is now 
exp[Xe(x2) — Xq^xi)], and the factor assigned to each turning point z is -^Fq{z), where 

F.(.) . (5.14) 
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Figure 9. Diagrammatic interpretation of (a) pr,e{x), and (b) pifi{x). Now the path 
is reflected at infinity. Each reflection at ±cxj gives a factor ±77. 



Since the expression (5.11) involves the state e '''"9%/2|0; 6'), it is convenient to define, 
in addition to (5.13), the quantities 

(0; e\aeaeU{x, -oo)e-^'*s"e/2|o. 



{0;e\U{x,-oo)e-'^^yo/'^\0;e) 
(0; ^|e^"«"«/2f/(^oo^ x)44|0; 6) 



(5.15a) 



(5.156) 



(0;^|e^«fl'^«/2[/(oo,x)|0;^) 

They can be interpreted as refiection coefficients including additional scattering at 
infinity (see figure 9). 

The expressions in the rotated frame corresponding to (4.12 a) and (4.126) are 
obtained by adding the subscript 6 to the scattering coefficients, and making the 
replacements ikx — )■ Xe^x) and f{z) — t- Fe^z): 

Te{x, x') = [1 - Ri,e{x, x')Rr,0{x', -00)] exp Xe{x) - Xe{x') - / Fe{z)Rrfi{z, -00) dz , 

L Jx' 

(5.16a) 

r 

Te{x, x') = [1 - Ri,e{oo, x)Rr,e{x, x')] exp Xe{x) - Xe{x') + / Fe{z)Ri^e{oo, z) dz . 

L J x' 

(5.166) 

From the derivation of (4.12), and from the diagrammatic interpretation of the quantities 
Prfi and pifi shown in figure 9, it is obvious that these expressions still hold when 
Rrfi{z, —00) and Ri^e{oo,z) are replaced by Pr,e{z) and pifi{z), respectively: 



Tg{x, x') = [1 - Ri,eix, x')pr,0{x')] cxp Xe{x) - Xe{x') - / Fe{z)prfi{z) dz 

I Jx' 

T0{x, x') = [1 - pi^e{x)Rr,e{x, x')] exp Xe{x) - Xe{x') + / Fe{z)pi^e{z) dz 

I Jx' 



(5.17a) 
(5.176) 



Comparing (2.5) with (5.11), we find that the generalized forms of (3.6a) and (3.66) 
are obtained by making the following replacements: First, the scattering coefficients for 
the interval (x', x) are replaced by the quantities with subscript 6: 

t{x,x') ^ Tg{x,x'), Rr{x,x') ^ Rrfi{x,x'), Ri{x, x') Ri^g{x , x') . (5.18) 
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Second, the reflection coefficients for semi- infinite intervals are replaced not by Rr^g and 
Ri^e but by pr,e and pi^e: 

Rr{z, -oo) Pr,e{z), Ri{oo, z) pi,e{z), (5.19) 

where z is either x or x'. This is because the state e~^"'e"'e/'^\0;9) appears in (5.11) 
instead of |0; ^). Third, since the operator a + a' is replaced by the right-hand side of 
(5.12), the expression [1 -|- Ri{oo, x)][l + Rr{x, — oo)] in (3.6) needs to be replaced as 



where 



[1 + i?i(oo, x)\[l + Rr{x\ -oo)] Ae{x, x'), 

Ae{x, x') = cos I - sin | + (cos | + sin |j pi^g 

cos I sin I (cos | - sin |) Pr,e{x') 



(5.20) 



X 



Making these replacements in (3.6) leads to 

Ag{x, x') Tg{x, x') 



G{x,x') 
G{x,x') 



[1 - pi^o{x)pr,o{^)] [1 - RlA^, x')pr,o{x')] ' 

A0{x,x') Te{x,x') 



[1 - pifi{x)Rrfi{x, X')] [1 - pifi{x')prfi{x')\ ' 

Substituting (5.17a) and (5.17&) into (5.22a) and (5.22&), respectively, we have 
A0{x,x') 



G{x,x') = 
G{x,x') = 



1 - pifi{x)pr,e{x) 
Ae{x,x') 



exp 



Xe{x) - Xe{x') - f Fe{z)pr,e{z) dz 

Jx' 

Xeix) - Xe{x') + f Fe{z)pi,e{z) dz 



. ^ exp 
1 - Pi,e{x')Pr,e{x') 

where 

X0{x) = ikxcosO + ^V{x) sinO, Fd{x) = iksinO + f{x) cosd. 

Since equations (5.23) hold for any 6*, we may let 6* ^ — 6* in (5.236): 

A_0{x,x') 



G{x,x') 



"Tvexp 



X^o{x) - X_o{x') + r F^g{z)pi,_e{z) dz 

Jx' 



(5.21) 

(5.22a) 
(5.226) 

(5.23a) 
(5.236) 

(5.24) 
(5.25) 



1 - Pi-e{x')pr-e{x') 

A symmetric expression is obtained by taking the geometric mean of (5.23a) and (5.25): 

Aq{x,x')A_q{x,x') 



G{x, x') 



[1 - pifi{x)pr,e{x)][l - pi-g{x')pr-e{x')]J 

X exp 1^ [Xe(x) + X_g(x) - Xe(x') - X_,(x')]| 



X 



exp |-i J ^ [Fe{z)pr,e{z) - F_0{z)pi-e{z)] dz^ . 



(5.26) 



Thus, we have obtained expressions for G{x,x') in terms of pr^ and pi g. 

The quantities Pr^o and pi^g can be expressed in terms of Rr{x, — oo) and Ri{oo,x), 
as we shall now see. Using (5.8), (5.9) and (5.4), we write (5.15a) as 



|g-W2(cos fa - sin |at)(cos fa - sin f at)[/(x, -oo)|0) 
(0|e-W2t/(a;^ -oo 



(5.27) 
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and 

[a^a\ e"'"^"/'] = (ry + 2r/a^a - r)^aa)e-'^''''/^ , 
equation (5.27) can be modified to the form 

2x(0|e-''''"/^aaf/(a;, -oo)|0) 

where we have also used the definition (5.10). Using the diagrammatic interpretation, 
we can easily see that 

(0|a2"f/|0) (2n)! 
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(5.28) 
(5.29) 

(5.30) 



• on 

{0\U\0) ^ 2^n\ '■■ 
(See [13] for an explanation.) Therefore, 



m\o) 



Hence we obtain 
^0|e-W2aa[/|o) 



d_ (o|e-W2^|o^ 



(0|e-W2f/|o) 
Substituting this into (5.30) gives 

Rr{x, — OO 



dr) 



log(l + r]Rr 



Rr 



l + TjRr 



V 



1 + rjRr{x, — oo) 
In the same way, we have 

Ri{oo,x) + Tj 



(5.3i; 



(5.32) 



(5.33) 



(5.34a) 



1 — r]Ri{oo, X 

Substituting (5.24) and (5.34) into (5.26) we obtain, after some calculation, 
n ^ ( [l + vRr{xm+vRi{x')] \ [1 + Ri{x)][l + Rr{x')]e''^^--'^ 

^^'^ ^ 1^[1 + vMxrn + VRr{^')]) [1 - Rr{x)Ri{x)Y/^[l - Rr{x')Ri{x')]y^ 

RJz) , RAz) 



(5.346) 



X exp 



X exp 



1, n, /-^ . / Rr(z) RAz) 



2^ ' 



dz 



(5.35) 



l + r]Rr{z) l + riRi{z) 

where Rr{x) and Ri{x) stand for Rr{x,—oo) and Ri{oo,x), respectively. This is the 
generalized form of (4.15). Equation (5.35) holds for any 9, i.e., for any real number 
r] = tan 6'. We recover (4.15) by setting 6* = (77 = 0). 
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6. Expression with 9 — 7r/2 

With 6 = ±7r/2 [t] = ±1), the right-hand side of (5.35) takes a form that does not 
include the function f{x) explicitly. In particular, we have a very simple expression 
with 9 — +7r/2. Let us define 

S(x-h)- ^^('^'^'^) Rr{x,-oo;k) 

^ ~ 1 + Ri{oo,x;k) 1 + Rr{x,-oo;k)' ^'^ 

Setting = 1 in (5.35), and using (6.1), we can write 

1 r /"^ 1 

G(x, x'\ k) — —i= exp ik(x — x') — ik S(z:k)dz . (6.2) 

^[1-S{x;k)][l-S{x';k)] I ' Jx' ' ' \ 

Thus, the Green function is expressed in terms of the single function S defined by (6.1). 

This expression is valid even when there arc bound states. The right-hand side 
of (6.2) becomes infinite when S{x; k) = 1. Note that S{x] k) — 1 \s equivalent to 
Ri{oo,x]k)Rr{x,—oo]k) — 1, which is obviously the condition for resonance. (This 
condition does not depend on x; if S{x;k) — 1, then S{x';k) — 1, too.) This means 
that /c^ is an eigenvalue of the Schrodinger operator if S{x;k) — 1. Otherwise, (6.2) is 
always finite §. 

7. Relation with the WKB approximation 

It is interesting to think about the connection between (6.2) and the WKB method. In 
the leading order WKB approximation, a wave function satisfying (1.1) has the form 



il)[x) ~ J- — exp 



X 



i I p{z) dz 



where p{x) is the local wavelength for the Schrodinger equation defined by 



(7.1) 



p{x) = ^k^-Vs{x). (7.2) 

The function G[x,x') given by (6.2), with fixed x', is an exact solution of (1.1) for 
X > x'. Comparing (7.1) with (6.2), we can see that the quantity k{l — S) in the exact 
expression corresponds to the local wavelength p in the WKB approximation. 

Let us first see how (7.1) can be recovered from (6.2). If f{x) is constant, say 
f{x) = c, then Rr and Ri take the form [9] 



Rr{x,—oo;k) — —Ri{oo,x;k) — '■ . (7.3) 

Suppose that, for a non-constant f{x), the refiection coefficients can be approximated 
by the same form as (7.3): 



ik — i\ k"^ — f^ix) —ik + i\/k'^—P(x) 
Rr{x,-oc;k)^ ^^^^^ Ri(oc,x;k)^ ^^^^ (7.4) 

§ If /(+oo) and /(— oo) are both finite, or if Imfc > 0, then \S(x; k)\ is finite for any x. In other cases 
it may happen that = oo for some x, but this causes no problems. 
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Substituting (7.4) into (6.1) gives the approximation for S, 

Six; k)^l-'^^k^-f\x). (7.5) 

If we assume that f{x) varies slowly so that /c^ — ~ /c^ — — /', then 

k[l - S{x] k)] ~ yjk'^ - Vs{x) = p{x). (7.6) 

This reproduces the WKB approximation (7.1). 

Let us make a more detailed comparison by considering the higher-order corrections. 
Since the WKB expansion is essentially a high-energy expansion, it can be compared 
with the expansion in terms of 1/k. The asymptotic expansion of the reflection 
coefficients in powers of 1/fc was studied in [11]. Using the formulas derived there, 
we can express the corrections to (7.4) as a series in powers of 1/k. We have 



Rr{x, — oo) 



Ri{oo, x) = —i 



' 1 

.k-y/k'- f 



f" 



f 



{2ikY {2ikY 

f ^ f" 



, 5ff'-f" 



+ 



{2ikY 
5ff' - f" 



+ 



(7.7a) 



(7.76) 



{2ik)^ ' {2ikf ' {2ikY 

The condition for the validity of (7.7) as an asymptotic expansion is discussed in [11]. 
Substituting (7.7) into (6.1), we obtain 

S{x) = 1 - ^^^^ + ^ + ^ [2/'/' - - 2//" - /'"] + ■ ■ ■ . (7.8) 

(A formula is available for the coefficient of l/A;" in the expansion (7.8) for an arbitrary 
positive integer n.) By using 



(7.9) 



we can rewrite (7.8) in terms of Vs as 



S{x) = 1 



+ 



(7.10) 



k 8k^ 32k^ ^ ^^^svs 

In the WKB method, on the other hand, the wave function incorporating the higher- 
order corrections is written as [18] 

1 



exp 



'W{x 

W ^Wo + Wi + W2 + 



/X 
W( 



z) dz 



where 

Woix) 



p{x), Wi{x) 



p" 



m 

4p2 ' Sp3 



/\2 



+ 



W2{X) 



16 p4 



+ 



etc. 



(7.11) 
(7.12) 

(7.13) 



The WKB expansion (7.12) is an expansion in powers of the constant h'^ (which we have 
set to be unity) which multiphes the fist term on the left-hand side of (1.1). It is easy 
to see that Wi = 0(1/ k^), = 0(1/A;^), etc as A; oo. (The terms of W2 omitted 
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in (7.13) are 0{l/k'^).) So we can rearrange (7.12) into an expansion in powers of 1/k. 
We have 



Wi 



/\2 



+ 



(4) 



+ 0(1/A;^), etc. (7.14) 



F 16A;5 32A;5 ' " " ^ 32k^ 

It is obvious that (7.10) is equivalent to (7.12) with (7.14), where W = 1 — kS. 

Let us next see how the Bremmer series can be described in our formahsm. For 
this purpose, it is convenient to make use of the rotation introduced in section 5 with 
an imaginary angle 6. We define 

ifix) 



9{x) = arctan ■ 



k 



(7.15) 



If / is a constant, the rotation (5.3) with this angle ^ is a transformation to the frame 
of coordinates in which no scattering takes place [17]. Using this 9{x), the quantity on 
the right-hand sides of (7.4) can be written as 

k-^l/'-f^x) 



The expressions (7.4) are exact if 9{x) is an x-independent constant. The corrections to 
(7.4) can be expressed as a series in powers of 9'{x) — {d/dx)9{x): 



Rr{x, — oo) = tan — sec — ^ 



2 J- 



dz9'{z)e 



2iAi 



+ 1 r dzi r dz2 r dz3 9'{zi)9'{z2)9'{zs)e'^' + 

o J—oo J —oo J —oo 



(7.17a) 



^{x) o9(x) 
Ri{oo,x) = — tan— sec — — 



dz9'{z)e 



-2iAi 



1 fOO POO fOO 

+ -/ dzi dz2 dz^9\z^)9\z2)9\z^)e-'^^ + 

O Jx J Zl J Zl 



where we have defined 



q(x) = ^k^-f\x), 
Ai = q{w) dw, 

rx rZi nzi 

q{w) dw + / q{w) dw + / q{w) dw + / q{w) dw. 

J Z-i J Z'2 J Z^ 



A 



(7.176) 

(7.18) 

(7.19a) 

(7.196) 



(Since it is not the purpose of the present paper to discuss the approximation methods 
for the reflection coefficients, we omit the explanation here. Let us only mention that 
(7.17) can be derived from equations (10.2) and (10.3) of [17].) From (7.17) and (6.1) 
we obtain 

q{x) 1 



S{x) 



k 



cos \9{x) + sin | 



dz 9'{z) exp 



2i I q{w) dw 



cos \9{x) — sin \9{x) 



roo rz 

/ ci^^'(2;)exp 2i / q{w)dw 

Jx L Jx 



+ 



(7.20) 
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(It is possible to construct a formula for the term of an arbitrary order in the expansion 
(7.20).) Just like (7.10), we can rewrite (7.20) in terms of p as 

p{x) 



S(x) = 1 



P[x) 
k 



+ 



p{x) 
~2k 



dz 



2k 
p'{z 



' dz^ 

oo p{z) 



exp 



2i / p(w) dw 



p{z] 



exp 



2i I p{w) dw 



+ 



Substituting (7.21) into (6.2), expanding the right-hand side, 
integration by parts, we obtain 



(7.21) 

and carrying out 



G{x,x') 



1 



^Jp{x)p{x'\ 



: exp 



p{z) dz 



1 + 



1 



dz 



p\z) 
p{z) 



exp 



/x' 
p{w) dw 



1 r°° , p'(z) 

- dz , , exp 

2 Jx p{z) ^ 



2i 



I P{w) 

Jx 



dw 



+ 



(7.22) 



which corresponds to the ordinary Bremmer series. Whereas the Bremmer series is an 
expansion of the wave function iIj{x)^ the expression (7.20) or (7.21) gives a similar 
expansion for the quantity corresponding to W{x) of equation (7.11). 

As we have noted, the WKB approximation (7.1) is obtained by replacing the 
Fokker-Planck potential V{x) by a linear function at each point x. Another possible 
approximation is to replace V{x) by a quadratic function at each x. The reflection 
coefficients for quadratic potentials can be exactly obtained [11]. By substituting these 
exact expressions into (6.2) with (6.1), we obtain an approximation for the Green 
function. In some cases, this approximation can be better than the WKB approximation. 

The methods related to the WKB approximation we have seen above is just an 
example of using (6.2) for approximate evaluation. In making an approximation, in 
general, it is easier to deal with the reflection coefficients than the Green function itself. 
For each approximation method for the reflection coefficients, the expression (6.2) gives 
the corresponding approximation for the Green function. 



8. Conclusion 



In this paper, we have derived some exact expressions for the Green function. A general 
symmetric expression is given by (5.35). Reflecting the symmetry of the Fokker-Planck 
equation, this expression includes an arbitrary parameter t]. The simplest expression 
(6.2) is obtained by setting r) = 1. Analytic properties of the reflection coefficients can 
be studied relatively easily. By using the expressions derived here, we can investigate the 
properties of the Green function on the basis of the analysis of the reflection coefficients. 
In particular, (6.2) is useful for studying the high-energy behavior of the Green function. 
It also serves as a starting point for various approximation methods. 

The reflection coefficients Rr{x, — oo) and Ri{oo,x) that appear in our expressions 
have been deflned by using the Fokker-Planck equation. Of course, this is not the only 
possible way of deflning reflection coefficients for semi-infinite intervals. It is also possible 
to express the Green function in terms of reflection coefficients deflned in a different 
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way, without using the Fokker-Planck equation. However, the resulting expressions 
become more comphcated if we use a different (inequivalent) definition of the reflection 
coefficients. 

References 

[1] Landauer R 1951 Phys. Rev. 82 80 

[2] Kira M, Tittonen I, Lai W K and Stenholm S 1995 Phys. Rev. A 51 2826 
[3] Bremmer H 1949 Physica 15 593 

[4] Bellman R and Wing G M 1976 An Introduction to Invariant Imbedding (New York: Wiley) 

[5] Chandrasckar S 1960 Radiative Transfer (New York: Dover) 

[6] Heinrichs J 1986 Phys. Rev B 33 5261 

[7] Kumar N 1985 Phys. Rev. B 31 5513 

[8] Rammal R and Doucot B 1987 J. Physique 48 509 

[9] Miyazawa T 1989 Phys. Rev. A 39 1447 
[10] da Luz M G E, Heller E J and Cheng B K 1998 J. Phys. A: Math. Gen. 31 2975 
[11] Miyazawa T 2006 J. Phys. A: Math. Gen. 39 7015 
[12] Riskcn H 1984 The Fokker-Planck Equation (Berlin: Springer) 
[13] Miyazawa T 1995 J. Math. Phys. 36 5643 
[14] Miyazawa T 2000 J. Phys. A: Math. Gen. 33 191 

[15] Stokes G C 1883 Mathematical and Physical Papers 2 (London: Cambridge University Press) 
[16] Zinn- Justin J 1989 Quantum Field Theory and Critical Phenomena (Oxford: Oxford University 

Press) 

[17] Miyazawa T 1998 J. Math. Phys. 39 2035 

[18] Bender C M and Orszag S A 1978 Advanced Mathematical Methods for Scientists and Engineers 
(Auckland: McGraw-Hill) 



